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Abstract 

Since Grover's seminal work, quantum search has been studied in great 
detail. In the usual search problem, we have a collection of n items 
xi,...,x n and we would like to find i : Xi = 1. We consider a new 
variant of this problem in which evaluating Xi for different i may take 
different number of time steps. 

Let ti be the number of time steps required to evaluate x%. If the 
numbers ti are known in advance , we give an algorithm that solves the 
problem in 0{^/t\ + t\ + . . . + t%.) steps. This is optimal, as we also show 
a matching lower bound. The case, when U are not known in advance, can 
be solved with a polylogarithmic overhead. We also give an application 
of our new search algorithm to computing read-once functions. 



1 Introduction 

Grover's quantum search algorithm |10| is one of two most important quantum 
algorithms. It allows to search a collection of n items in 0{-^/n) quantum steps. 
This gives a quadratic speedup over the exhaustive search for a variety of search 
problems [2]. 

An implicit assumption is that examining any two items can be examined 
in the same number of time steps. This is not necessarily true when Grover's 
algorithm is applied to a specific search problem. It might be the case that some 
possible solutions to the search problem can be checked faster than others. 

Let ti be the number of time steps required to check the i th solution. Clas- 
sically, searching for an item i : Xi = 1 requires time 0(ii + . . . + t n ). A naive 
application of Grover's search would be to use 0(y/n) steps, with the maximum 
possible query time t max = max^ ti in each step. This gives a 0{\fnt max ) time 
quantum algorithm. 

In this paper, we give a better quantum algorithm. We consider two settings: 

1. The times ti are known in advance and can be used to design the algorithm; 

2. The times ti are not known in advance. The algorithm learns ti only if it 
runs the computation for checking the i th item for ti (or more) steps. 
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For the first setting, we give a quantum algorithm that searches in time in 
time O(VT) where T = t\ + ... + t\. For the second, more general setting, 
we give an 0(\/T log 2 TTog 2 logT) time quantum algorithm. We show a lower 
bound of f2(v / T) for the first and, hence, also the second setting. 

We give an application of our search algorithm, to computing read-once 
Boolean functions. A Boolean function f(x±, . . . , Xjv) is read-once if / has a 
formula (consisting of AND, OR and NOT operations) in which every of the 
variables x\, . . . , xn appears at most once. We show that any read-once Boolean 
function for which the depth of the read-once formula is d can be computed 
using 0(VN "log^ 1 N) queries. Previously, such algorithm was only known for 
the case of balanced AND-OR trees PITT]. 

The model in which queries to different x% take different time has been 
previously studied by H0yer et al. ^2] wn ° proved composition theorems for 
quantum lower bounds in a similar model. Our paper appears to be the first to 
study the complexity of quantum search in such model. 

2 Model 

We would like to model the situation when the variable x% is computed by an 
algorithm Ai which is initialized in the state |0) and, after U steps, outputs the 
final state \xi)\ipi) for some unknown \ipi). (For simplicity, we assume that Ai 
always outputs the correct Xi.) In the first £j — 1 steps, Ai can be in arbitrary 
intermediate states. 

Our goal is to find i : Xj = 1. (We sometimes refer to i : Xi = 1 as marked 
items and i : Xj = as unmarked.) Our procedure A can run the algorithms Ai, 
for some number of steps t, with Ai outputting X, if ij < t or "the computation 
is not complete" if t% > t. The computational cost is the amount of time that is 
spent running algorithms Ai- Any transformations that do not involve Ai are 
free. This is a generalization of the usual quantum query model. 

For completeness, we include a more formal definition of our model in the 
appendix. Our algorithms, however, can be understood with just the informal 
description in the previous two paragraphs. 

Known vs. unknown times. We consider two variants of this model. In 
the "known times" model, the times t±, . . . ,t n are known in advance and can 
be used to design the algorithm. In the "unknown times" model, t\, . . . , t n are 
unknown to the designer of the algorithm. 

3 Methods and subroutines 
3.1 Amplitude amplification 

Amplitude amplification is a generalization of Grovcr's quantum search al- 
gorithm. Let 

sma|l)|'!/>i) + cosa|0)|^>o) (1) 
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be the final state of a quantum algorithm A that outputs f with probability 
sin 2 a — 8. We would like to increase the probability of the algorithm outputting 
1. Brassard et al. [7] showed that, by repeating A and A" 1 2m + 1 times, it is 
possible to generate the final state 

sin(2m+ l)a|l)|V>i) + cos(2to + l)a|G)|V>o)- (2) 

In particular, taking m — 0(^=) achieves a constant probability of answer 1. 

We use a result by Aaronson and Ambainis who gave a tighter analysis 
of the same algorithm: 

Lemma 1 [|J? Let A be a quantum algorithm that outputs a correct answer and 
a witness with probability 1 8 < e where e is known. Furthermore, let 

to < = . (3) 

4 arcsin 2 

Then, there is an algorithm A' which uses 2m + I calls to A and A^ 1 and 
outputs a correct answer and a witness with probability 

S new >(l- {2m + 1)2 s) (2m + 1) 2 <5. (4) 



The distinction between this lemma and the standard amplitude amplifica- 
tion is as follows. The standard amplitude amplification increases the probabil- 
ity from 8 to in 2to + 1 = 0(-^=) repetitions. In other words, 2to + 1 rep- 
etitions increase the success probability Sl((2m + l) 2 ) times. Lemma ^ achieves 
an increase of almost (2to+ l) 2 times, without the big- SI factor. This is useful if 
we have an algorithm with k levels of amplitude amplification nested one inside 
another. Then, with the usual amplitude amplification, a big-Si constant of c 
would result in a c k factor in the running time. Using Lemma avoids that. 

We also need another fact about amplitude amplification. 

Claim 1 Let 8 and 8' be such that 8 < e and 8' < e and let to satisfy the con- 
straint (0). Let p(8) be the success probability obtained by applying the procedure 
of Lemma^to an algorithm with success probability S. Lf 5' < 5 < cd' for c> 1, 
then p(5') <p{5) < cp{8'). 

Proof: Because of equations (Q, 

p(8) = sin 2 ((2TO + 1) arcsinv^). 

Let 7 = arcsin y/d and 7' = arcsin y/S'. Then, we have to prove that sin 2 7' < 
sin 2 7 < csin 2 7' implies sui 2 (2to + 1)7' < sin 2 (2TO + 1)7 < csui 2 (2to + 1)7'. 

Because of 8 < e and 8' < e, we have yj~8 < yfe and y/W < y/e. Together with 
J2J, that means that (2m + 1) arcsin y/S < ? and (2to+ 1) arcsin y/W < ?. Since 

1 HI requires the probability to be exactly e but the proof works without changes if the 
probability is less than the given e. 
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sin is an increasing function on [0, sin 7' < sin 7 implies sin (2m + 1)7' < 
sin 2 (2m + 1)7. 

To prove the other inequality, we consider the function f(x) = ■ 
It suffices to show that /(x) is non- increasing on [1,2m + 1]. (That implies 

snr((2m4-l)7 ) — sin y ' 

cos(x7) sin(x7 / ) — sin(x7) cos(x7') sin #(7' — 7) 



/'(*) = 



sin (xY) sin (27') 



This is non-positive, as long as x(j' — 7) € [^-,0]. This is true, because 
^(Y — 7) > — ^7 > —(2m + 1)7 > — ^ and 2(7' — 7) < follows from 7' < 7 
(which follows from 5' < 8), | 

3.2 Amplitude estimation 

The second result that we use is a version of quantum amplitude estimation. 

Theorem 1 ^ There is a procedure Est-Amp(yl, M) which, given a quantum 
algorithm A and a number M , outputs an estimate e of the probability e that A 
outputs 1 and, with probability at least we have 



I -I ^ \/max(e(l - e),e(l - e)) tt 2 

The algorithm uses M evaluations of A. 

We are interested in a slightly different type of error bound. We would like 
to have |e — e| < ce for some small c > 0. 

Theorem 2 There is a procedure Estimate(_4, c,p, k) which, given a constant 
c, < c < 1 and a quantum algorithm A (with the promise that the probability 
e that the algorithm A outputs 1 is either or at least a given value p ) outputs 
an estimate e of the probability e such that, with probability at least 1 — ^f, we 
have 

(i) |e - e| < ce if e > p: 
(li) e = ife = 0. 
The procedure Estimate(_4, c,p, k) uses the expected number of 



6 I k[ 1 + loglog- 



p J V max(e,p) 



evaluations of A. 
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1. Let M = 2; 

2. Repeat: 

(a) Let e be the estimate output by repeated Est-Amp(„4, M). 



(b) If 271- ^^ — + j-p < tie, stop and output e as the estimate. 

(c) M = 2 * M. 



until M > M moa! where M mnx = 



8tt 



Algorithm 1: Procedure Estimate 



Proof: We can increase the success probability of Est-Amp(_4, M) to at least 
1 - WT^Tm (where M max = -M ), by repeating the algorithm t = 

& lOgJIJmax CyJ (1 — C)p 

0((1 + log log ~)fe) times and taking the median of the results. 

The procedure Estimate calls the repeated Est- Amp at most \ogM max 
times. Since each call of Est- Amp produces the correct answer with probability 

at least 1 — 2k log \j , the probability that all calls to Est- Amp produce correct 

results is at least 1 — In this case, Estimate is always correct, because by 

Theorem ^ the error |e— e| is at most + -™ and Estimate only stops 

when this quantity becomes less than ce. It remains to bound the number of 
times Estimate calls A. 
If M > 4ff , then, 

<V(l-c)e 



2 J^E± + 2 4 < ?£Z± + < C (i - c)e . (5 ) 

M Af 2 - 2 16 ~ v ; w 

Then, e > (1 — c(l — c))e > (1 — c)e. Therefore, the quantity of equation JSJ is 
less than or equal to ce. Hence, if M > — 4?r , then the condition in step !2bl 

is satisfied and the algorithm stops. Since M is doubled in every iteration, the 

final value of M is M < — , 87r The algorithm A is repeated 

<V(i-c)« 

M Q t M t „ „ 16tt 

M i + + + . . . < 2M t < t 
2 4 <V(l-c)e 

times. 

If e > p, the algorithm must stop with M being at most — jSS, . If that 

ci/(l-c)p 

does not happen, we can conclude that e = 0. The number of repetitions of A 

in this case is at most — ^22 — t. I 
<V(i-<0i> 
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4 Search algorithm: known running times 



Theorem 3 A collection of n items with times t\,...,t n can be searched in 
time 



0\Jtj + ti + ... + t; 

Proof: The basic idea is to subdivide the items into groups so that all items in 
one group have similar times ti (e.g. < t\ < t max for some t max ). We can 

perform the standard Grover search in a group in time s = 0(^/lt max ) where / 
is the size of the group. We then observe that 



s 2 = 0(lt 2 max ) 



°(? tf ) 



with the summation over all items i in the same group. By summing over all 
groups, we get 



where j on the left ranges over all groups. Let k be the number of the groups 
that we have. If we have a search algorithm that searches k items in time 



k j 



we can then substitute the algorithms for searching the k groups instead of the 
k items and obtain a search algorithm for n items that runs in time 



o(y/tl + ... + t; 



We then design a search algorithm for k items in a similar way. 

The simplest implementation of this strategy gives an algorithm with log* n 
levels of recursion and running time 



o c lo «*V*? + *! + ••• + *« 



due to the reduction from n items to k items losing a constant factor every 
time it is used. The c log ™ factor can be avoided, by a more sophisticated 
implementation of the same idea, which we describe below. 

We first restrict to the case when there is exactly one marked item. The 
general case can be reduced to this case with a constant factor overhead, by 
running the algorithm on all n elements, a random set of ^, a random set of j, 
etc. As shown in pQ, there is a constant probability that at least one of those 
sets contains exactly one marked item. The expected running time increases by 
at most a constant factor, because of the following lemma. 



G 



Lemma 2 Let S be a uniformly random set of elements of {1,2, ... ,n} . 
Then, 



E 




< 



2i/ 2 • 



E n 

«£{!,.. .,n} 



Proof: By concavity of the square root function, 



E 




E 



ies 




Therefore, the reduction from the general case to one marked item case 
increases the bound on the number of queries by a factor of at most 



1 + 



1 

21/2 



. . . < 



Second, we introduce a generalization of the problem in which the algorithm 
Ai for the marked i returns the correct answer with a probability at least pi, 
instead of a certainty. More formally, 

• if Xi — 0, the final state of the algorithm Ai is of the form |0}|V>o)- 

• if Xi = 1, the final state of the algorithm Ai is of the form a|l)|V>i) + 
Vl — a 2 |0)|^o)j where pi < \a\ 2 < d • Pi, for some constant d > 1. 

The probabilities p\,...,p n and the constant d are known to us when we design 
the algorithm, just as the times t\, . . . ,t n . (Knowing both the success proba- 
bility and the running time may look quite artificial. However, we only use the 
"known success probability" model to design an algorithm for the case when all 
Ai return the correct answer with certainty.) 

We claim that, in this case, we can search in time 




Our main theorem now follows as the particular case p\ 
main part of our proof is 



p n = 1. The 



Lemma 3 There exists k — O (log 3 n log log n) with the following property. As- 
sume that there is a search algorithm for k items with some fixed d > 1 that 
works in time at most 



I 2 . 

92 



. + 
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for any given times s%, . . . ,Sk and probabilities q%, . . . , q^. Then, there exists a 
search algorithm for n items with d! = (l — O ( ^ - \ d instead of d that works 
in time at most 

1 



C [1 + 




Jogn 

for any given times t±, . . . , t n and probabilities pi, . . . ,p 



Proof: Let Tn be the maximum of -43= . . . . We first check all items 
with < - n sequentially. To check item i, we just amplify the success 
probability of Ai to by the standard amplitude amplification, in 0(^/p7) 
steps. Therefore, the number of steps for checking item i is 



Pi J \VPiJ \nlogn 

The time for checking all such items is at most the number of such items 
times — — which is of the order at most 

n log n 




log n log n 



Next, if ^ < we choose m so that < (2m+l) 2 Pi < (Such 

choice of m always exists, because, if (2m + l) 2 pi < § t * - , then 

ox2 (2m + 3) 2 1 „ 1 1 
(2m + i) 2 Pl < )- -4 7T1 ^ 9 7n ^ 



(2m+l) 2 91ogn 91ogrt logn 

Therefore, it suffices to choose the smallest m for which (2m + l) 2 pi > 9l ^ gn ■) 
We then apply Lemma^ If the success probability is pi, it increases the success 
probability to p'^ while increasing the running time (2m + 1) times. By Claim 
if the success probability is between P i and d- pt, it increases to a probability 

m t 2 

between and d ■ p^. By Lemma ^ the ratio increases at most 



3 

times. After that, we have 



j _ (2m+iP " ! _ 1_ yi ogn 
^ o Vt 3 log n \ o 



Claim 2 (a) pa < Pi < I, where po — 1 — 



91 



1 / n log n — L — u 
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Proof: (a) Let p\ be the value of pi before the amplification. If p\ > 
then Pi = p\. If p\ < then 

p , = f 1 _2!!i±iMV 2m+im >f 1 - 1 A 1 



3 / \ 3 log nj 9 log n 

(b) The first inequality follows from -^fe= > - - and pi > po- The second 
inequality follows from To = max -Ji= and Pi < 1. | 

We partition the intervals [^^, T ] and [p , 1] into subintervals of the form 
[T',T"] and [p',p"] with T" < (1 + and p" < (1 + It suffices to 

have 0(log 2 n) intervals for time and O(lognloglogn) intervals for the success 
probability. The overall number of intervals is m = 0(log 3 n log log n). 

We define new search algorithms A'- , where j ranges over the pairs of intervals 
[T', T"] and \p',p"]. Let [Tj, Tj'] and fy'pp"] be the time and probability intervals 
for A'j and let Sj be the set of i such that Tj < U < T'j and p' 3 < p t < p" . The 
algorithm A'j picks i G Sj uniformly at random and then runs Ai . 

Let Sj denote the running time of A'j. Then, Sj < T". The success proba- 
bility of A'j is in the interval [qj, d(l + y^)qj\ where qj = t^t. We now relate 
Sj and qj to Tj and pj : 

yji>is,i<H2! = ('i + o (+)) >(i + o(')) M. 

p 3 \ \lognJJ lJl p'j \ \\ognJJ qj 

By summing over all pairs of intervals j, 

s 2 s 2 / / 1 \\ /t 2 t 2 



9i <?fc V J \Pi P 

We now apply the search algorithm for fc items to A'i, . . . , A' k . § 

To obtain Theorem we repeatedly apply Lemma until the number of 
items becomes less than some constant no- That happens after 0(log* n) appli- 
cations of Lemma |3] 

Let ti,...,t n and pi, . . ., p n be the times and probabilities for the final 

n < n-o items. After that, we just amplify the success probability of every item 

t 2 

to r2(l) (which increases each ^- by at most a constant factor, as discussed in 
the proof of Lemma EJ. We then search n items in time 0(^/nniaxi U), using 
the amplitude amplification, with max^ £j steps for evaluating any of the items 
i. Since pi = fi(l) and n < no where no is a constant, we have 




\/nmaxti — O(maxii) = O yy t\ + . . . + t 

0(log* n) applications of Lemma |3] increase the time by a factor of at most 
( 1 + 0(^)) l ° s " n = l + o(l). I 
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5 Application: read-once functions 



A Boolean function f(xi, . . . ,%n) that depends on all variables xi, . . . ,Xn is 
read-once if it has a Boolean formula (consisting of ANDs, ORs and NOTs) 
in which every variable appears exactly once. A read-once function can be 
represented by a tree in which every leaf contains Xi or NOT Xi and every 
internal vertex contains AND or OR. 

Barnum and Saks 0] have shown that, for any read-once /, Q(y/~N) queries 
are necessary to compute / in the quantum query model. This bound is known 
to be tight for a special class of read-once functions: balanced AND-OR trees. 
A balanced AND-OR tree is a read once function represented by a depth-d tree 
in which each internal node has \/~N children. Nodes on the even levels are AND 
nodes, nodes on the odd levels are OR nodes (or opposite). Hoyer, Mosca and 
de Wolf jllj have shown that, for any constant d, the function corresponding 
to the AND-OR tree of depth d can be evaluated with 0(y/~N) queries. This 
improved over an earlier 0{VN log 4 ' 1 N) query algorithm by Buhrman, Cleve 
and Wigderson (§] . Both of those algorithms depend on the fact that every node 
on the same level of the tree has an equal number of children. 

We give the first quantum algorithm for the general case, when the number 
of children may vary for different nodes at the same depth. 

Theorem 4 Any read-once function f(xi, . . . , xn) of depth d can be computed 
by a quantum algorithm that uses 0(VN log*- 1 N) queries. 

Proof: By induction. The base case, d = 1 is just the OR (or AND) function 
which can be computed with 0(y/N) queries using Grover's search to search for 
i : Xi = 1 (or i : Xi = 0). 

For the inductive case, assume that / is represented by a depth-d tree with 
OR at the root. (The case when the root contains AND is similar.) Let n be 
the number of vertices on the level 1 (that is, the number of children of the root 
vertex) and tj be the number of vertices in the subtree rooted in the i th level- 1 
vertex. By re-ordering the variables, we can assume that 

f(xi, . . .,x N ) = V" =1 / i (a; tl+ ... + t j _ 1 +i, ■ • .,x tl+ ... +f4 ). 

To compute /, we have to determine if there exists i £ {1, . . . ,n} for which 
fi = 1. By the inductive assumption, there is an algorithm that computes fi 
using 0(-\/U log d ~ 2 U) — 0(y/U\og d ~ 2 N) queries. We repeat this algorithm 
O(logiV) times to increase the probability of correct answer to at least 1 — -^j. 
Let Ai be the resulting algorithm and 1$ = 0(v^7log rf_1 N) be the number of 
queries in A4. 

We now apply Theorem [3] to A\, . . ., A n . This gives an algorithm which 
uses 

O Ut* + T% + ... + t£\ = 0{\og d ~ l NVt!+t 2 + ...+t n ) = OiVNlog^ 1 N) 
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queries. Since we are applying Theorem to Ai, . . ., A n which are incorrect 
with a small probability, we have to bound the error probability for the resulting 
algorithm. 

Let A[, . . ., A' n be the "ideal versions" of Ax, . . ., A n . If the final state of 
At is 

a\a)\il> a ) + Vl - a 2 \l - a)|^>i_ a ), (6) 

where a is the correct answer (the value of fi), then the final state of A\ is 
| a) IV'o)- (A'i can be obtained by composing Ai with a transformation that maps 
the state © to |a)|^ a ).) 

Given the "ideal algorithms" A' i: the algorithm of Theorem [3] would output 
the correct answer with a constant probability (e.g., at least 2/3). Since each 
Ai outputs the correct answer with probability at least 1 — -J^ , replacing Ai by 
A^ in one time step changes the state of the algorithm by at most O(^) (in the 
I2 norm). Replacing Ai by A\ in every time step changes the state by at most 

in I2 norm. Therefore, the success probability will still be | — o(l), even if the 
actual Ai , . . . , A n are used. § 



6 Search algorithm: unknown running times 

In some applications, it may be the case that the times ti are not known in 
advance. We can also solve this case, with a polylogarithmic overhead. 

Theorem 5 Let e > 0. There is an algorithm that searches collection of n 
items with unknown times ti,...,t n and, with probability at least 1 — e, stops 
after 

O (T log 2 T log 2 logT) 
steps, where T = yjt\ + t\ + . . . + t\ . 

Proof: Again, we assume that there is exactly one marked item. (The reduction 
from the general case to the one marked item case is similar to one in the proof 
of Theorem El) 

Let St be the set of items such that Xi = 1 or ti > 2* and let rit = \St\- Our 
main procedure, algorithm |21 defines a sequence of algorithms B\, . . ., Bi. The 
algorithm Bj, with some success probability, outputs a bit 1 and, conditional 
on output bit 1, it also outputs a uniformly random index i G Sj. To avoid 
the problem with accumulating constant factors (described after Lemma U, we 
make the success probability of Bj slightly less than 1. 

Lemma 4 Assume that the constant D in stevs Uial and \2ct satisfies D < -^=- 
Then, with probability 1 — e, the following conditions are satisfied: 
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1. Set j = 1. Define B\ as the algorithm that just outputs 1 and a uniformly 
random i G {1, . . . , n}. 

2. Repeat: 

(a) Use the algorithm Bj to generate k = 2log(D(j + 1)) samples 
ii, . . . ,ik of uniformly random elements i € Sj. Run 2 J+1 steps of the 
query procedure on each of ix, ■ ■ ■ , ih- If %i — 1 for one of samples, 
output i and stop. 

(b) Let Bj +1 be an algorithm that runs Bj once and, if the output bit 
is 1, takes the output index i and runs 2- 5 " 1 " 1 steps of the checking 
procedure on i. If the result is Xi — 0, outputs 0. Otherwise, it 
outputs 1 and the same index i. 

(c) Let p = Estimate(^ +1 , c, -±-,2 log(D(j + 1))). If p = 0, output "no 
i : x % = 0" . 

( d ) If P> 9hb> let ^+i heB 'j+v 

(e) If p < 9 t p , let Sj+i be the algorithm obtained by amplifying 

2m + 1 times, where m is the smallest number for which ttt^ — < 

1 9 log n — 

(2m + l) 2 p < • (Such choice of m always exists, as described in 
the proof of Lemma [3]) 

(f) Let j=j + l. 

Algorithm 2: Search algorithm for unknown t\, . . . , t n 



(a) Estimates p are accurate within an multiplicative factor of (1 + c); 

(b) If Bj is defined, then U > 2 j l for at least values i € {1, . . . , n}. 

Proof: (a) The probability of error for Estimate is at most jw^qrrp ■ By 
summing over all j, the probability of error for some j is at most 

1 1 1 7T 

i— 1 

which can be made less than -| by choosing D < -?=. 

(b) By definition, jS^-i is the set of all i with the property that either = 1 
or ^ > 2 J_1 . Let 5 be the set of i with x. t = 1 and U < 2 3 1 . If |5| < |nj_i, 
(c) is true. Otherwise, the probability that each ij generated in step |2a| does 
not belong to S is less than ^ . If one of them belongs to S, algorithm [21 stops 
without defining Bj. The probability that this does not happen (i.e., all ij do 
not belong to S) is less than (|) fe = D -i(J +1 yi ■ We can make this probability 
arbitrarily small similarly to part (a). | 
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We now bound the running time of algorithm [3 under the asumption that 
both conditions of Lemma0]are true. For that, we first bound the running time 
of the algorithms Bj and then the total running time of algorithmic We assume 
that both conditions of Lemma 01 are true. 

Let pj be the success probability of Bj and p'j be the success probability of B'-. 
Let ru,i be the number of times step|2e]is performed, for j S {k, k+ 1, . . . , I— 1}. 



Lemma 5 The running time of Bj is at most 



lt n'''./M + ff 1 + Ar.Sd 



log n/ y rij logn/ V Pj>-irij 



V (7) 



for some constant C. 



Proof: By induction. The base case is easy. Since pj — 1 and rij — n, the 
expression (Q is just equal to 1, which is also the running time of B\. 

For the inductive case, we first consider the running time of B'j +1 . It can be 
decomposed into two parts: the running time of Bj and the running time of the 
2 J+1 -step checking procedure. The running time of Bj is described by equation 
101 . We have p'j +1 = Pj " J+1 . Therefore, we can rewrite J7J) as 



logn/ W nj+i V l°g» 1 / f P]'-i« 3 +i 

The time for the checking procedure is just 2 3+1 which is equal to P,+1 " J 2 J+1 
(since = 1). Therefore, the running time of 6' , n is 



+ MtE+Y, (i + °)" J J^v. (8 ) 

logn/ y n J+ i ^ V log"-/ \pj>-in j+ i 

If step !2dl is performed, then Bj+i = B'j +1 , Pj+i = p'j+i, Tf,j = r j',j+i and the 
expression (JHJ is the same as {7} with j + 1 instead of j. 

If the step[2e|is performed, the running time of Bj+i is (2rn + 1) times the 
running time of B'j +1 . The success probability is 



ft+1 >( 1 _(^ p ; +1 ) (2m + lW+I >( 1 __i n 



(2m + 1) V 



Therefore, 



2m+l<fl + -^V/^±i (9) 
V lognJ^p' j+1 

for some constant C . Multiplying (JSJ) by 2m + 1 and applying 10 completes the 
induction step. | 
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Lemma 6 For all rjj> = 0(\ogn). 

Proof: We consider the ratio qj — We have 5i = ^ and q,j < 1 for all j 
(since pj < 1 and rij > 1). 

Next, we relate qj and qj+%. We have = If step [23| is applied, 

Pj+i = p'j+i and qj+i = = qj. If stepEilis applied, 



Therefore, qj+i > 9(1 — 3 i 1 g „ )ffj- This means that qy > (9 — 3 log n) r i.i' qj. 
Together with qji < 1 and qj > q\ > i, this implies = O(logn). | 
The expression of Lemma can be upper-bounded by 

Lemma 7 T/ie running time of Bj is 



Proof: We look at each of the components of the sum J7J) separately. Consider 
a term 

\ + -±-Y\m^2f. ( io) 

logn/ y Pj'-iUj 

Because of Lemma El the first multiplier is bounded from above by a con- 
stant. Since Pj'-i > g7ogri (similarly to Claim |2J), we can upperbound 11011 by 

O^logn "^ 1 2 3 ). Let k be the number of i S {1, . . . , n} for which £j > 2 J _1 . 

By Lemma 01 fc > 3 1 -1 and 

+ if + . . . + ^ > /c2^' ^ > 
This means that each term in <|10[) is at most 



o 



H\ + q + ...+ti 



The lemma follows by summing over all j terms in (01 . | 

We now bound the overall running time. To generate a sample from Sj, one 
needs O(yTogn) invocations of Bj (because the success probability of Bj is of 
the order ^(j^-^))- Therefore, we need O(^lognlogj) invocations to generate 
O(logj) samples in step|2a| By Lemma [7| that can be done in time 



O jlogjlogm 



iti + 4 + ... + ti 
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For each of those samples, we run the checking procedure with 2 J+1 steps. That 
takes at most twice the time required by Bj (because Bj includes the checking 
procedure with 2- 7 steps). Therefore, the time for the 2 J+1 checking procedure 
is of the same order or less than the time to generate the samples. 

Second, the success probability estimated in the last step is of order Pj " J+1 = 
^( nOogn ) - By Lemma |5J it can be estimated with 



O log j log log Ha 



' rij log n 
n J+1 



invocations of Bj, each of which runs in time described by Lemma 

Thus, the overall number of steps in one loop of algorithm is of order at 
most 

'jlogjlogn j log j log n log log n x 



t\ + t\ + . . . + 1\ 



Since rij > 1 and nj+i > 1, this is of order 



O I sjt\ + t\ + ■ . ■ + tl j log j log n log log n 

be the maximum of ti, . . ., t n . Then, the maximum value of j is at 
most \\og(t max + 1)]. Therefore, the number of steps used by the algorithm[3 
is 



[ V 1 1 + 2 + ■ ■ ' + log n log log 71 log tmax log log tm ' 

The theorem now follows from n < \[T and t max < VT, where T = t\ + 1\ 
... + tl | 



7 Search lower bound 

Theorem 6 For any positive integers ti, . . . ,t n> searching a collection of n 
items that can be checked in times t\, .. . ,t n requires time Cy/t\ +t 2 , + ... + t 2 l , 
for some constant c > 0. 

Proof: Let t\ be the maximum integer such that [f \f^\ + 1 < U (with ^ = 1 
if the maximum integer is 0). We consider searching m = t[ + . . . + t' n elements 
xi, . . . ,x m g {0, 1} in the standard model (where every query takes 1 step), 
with the promise that there is either or 1 element j ' : Xj — 1 . By lower bound 
on quantum search, c'-^m queries are required to distinguish between the case 
when there are elements j : Xj = 1 and the case when there is 1 element 
j : Xj = 1, for some constant c'. 

We subdivide the inputs xi, . . . , x m into n groups 5*1, . . ., S n , with ^, . . . , t' n 
elements, respectively. Let y.i = 1 if there exists j € Si with Xj = 1. Since there 
is either or 1 element j : Xj — 1, we know that there is either or 1 element 
i : yi = 1. We will show 
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Lemma 8 There is an algorithm that implements the transformation \i) — > 
/ or some states \ipi), using U queries. 

Let A be a search algorithm for search among n items that require times 
ti, . . . ,t n and let t' be the number of steps used by A. Then, we can substitute 
the algorithm of Lemma [S] instead of the queries yi. Then, we obtain an algo- 
rithm A' that, given X\, . . . , x n , asks t' queries and distinguishes whether there 
is exactly 1 item i : = 1 (and, hence, 1 item j : Xj — 1) or there is no items 
i : j/j = (and, hence, no items j : Xj = 1). Hence, 

t! > c'Vn = cV*i + ---+^- 
We now bound t\ in terms of tj. By definition of f^, we have 

+ 1< 

This means that t' t > $(t t - 2) 2 . If U > 3, then t; - 2 > ^ and t- > ^t?. If 
t ; < 3, then ^ > 1 > ^tf. Therefore, 

. H[6 4c' I 

t' > dVti + -• + *'„ > c 'y 9^2 (*? + ' ' • + f ™) = 3^V*i + ' • ■ + 

This means that the theorem is true, with c = |p. It remains to prove Lemma 

El 

Proof: [of Lemma [H] To simplify the notation, we assume that the group Si 
consists of variables xx,...,Xt i . If U = 1, then yi = x\ and we can just query 
x\. This produces the required transformation \i) — > \i)\yi). 

For the U > 1 case, we have to search ti items x\, . . . ,x ti for an item j : 
Xj = 1, if we are promised that there is either or 1 such item. There is a 
modification of Grover's algorithm which succeeds with probability 1, using at 
most [~f\/i7| queries jZj. 

The result of Grover's algorithm is: 

• the state where j is the index for which Xj — 1, if such j exists; 

• the superposition -7= X^=i U)j otherwise. 

With one more query (which queries the index j), we can determine the value 
yi = Xj (which is 1 in the first case and in the second case). | 

8 Conclusion 

In this paper, we gave a quantum algorithm for the generalization of Grover's 
search in which checking different items requires different times. Our algorithm 
is optimal for the case when times U are known in advance and nearly optimal 
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(within a polylogarithmic factor) for the general case. We also gave an appli- 
cation of our algorithm to computing read-once Boolean functions. It is likely 
that our algorithms will find other applications. 

While we have mostly resolved the complexity of search in this setting, the 
complexity of other problems have not been studied at all. Of particular inter- 
est are problems which are frequently used as a subroutines in other quantum 
algorithms (for such problems, there is a higher chance that the variable-time 
query version will be useful). Besides the usual quantum search, the two most 
common quantum subroutines are quantum counting [5] and A:-item search (a 
version of search in which one has to find k different i for which Xi = 1). Ele- 
ment distinctness [21 E| has also been used as a subroutine, to design quantum 
algorithms for the triangle problem |14j and verifying matrix identities jBJ I13j . 

Acknowledgments. I would like to thank Robert Spalek and Ronald de 
Wolf for the discussion that lead to this paper. 
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A Formal definition of our model 

To define our model formally, let be the j th step of Ai. Then, 

A- - A {u) A {ti ~ 1] A {1) 

We define ^4-*^ = I for t > ti. We regard the state space of Ai as consisting of 
two registers, one of which stores the answer (c £ {0, 1, 2}, with 2 representing 
a computation that has not been completed) and the other register, x, stores 
any other information. 

The state space of a search algorithm is spanned by basis states of the 
form \i, t, t r , c, x, z) where i £ {!,..., n}, t,t r £ {0,1, ... ,T} (with T being the 
number of the query steps in the algorithm), c £ {0, 1,2} and x and z range 
over arbitrary finite sets, i represents the index being queried, t represents the 
number of the time step in which the query for Xi started and t r is the number of 
time steps for which A will run the query algorithm Ai. c is the output register 
of Ai and x holds intermediate data of Ai. Both of those registers should be 
initialized to |0) at the beginning of every computation of a new X{. z contains 
any data that is not a part of the current query. 

We define a quantum query algorithm A as a tuple (Uq, . . . , Ut) of uni- 
tary transformations that do not depend on x\, . . . ,x n . The actual sequence of 
transformations that is applied is 

Uo, Qi, Ui,Q2, ■ ■ ■ , Ut-i, Qt, Ut, 

where Qj are queries which are defined below. This sequence of transforma- 
tions is applied to a fixed starting state \ip start), which consists of basis states 
\i,0,0,c,x,z). 

Queries Qj are defined in a following way. If j < t + t r , we apply A± *' to 
|c) and \x) registers. Otherwise, we apply /. We call the resulting sequence of 
queries Qi, Q2, ■ ■ ■ generated by transformations A\. We call Q\, Q2 a valid 
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sequence of queries corresponding to x\ , . . . , x n if it is generated by A\ satisfying 
the following constraints: 

1. For t < U, AjA*- 1 . . . Aj\0) is of the form |2}|V>) for some 

2. For t = U, A\A l r X ■ ■ ■ Aj\0) is of the form \xi)\ip) for some \ip). 

Uj can be arbitrary transformations that do not depend on xi, . . . , x n . 

An algorithm (Uq, • ■ • , Ut) with the starting state start) computes a func- 
tion f(x\, . . . , x n ) if, for every x\,...,x n <G {0,1} and every valid query se- 
quence Qi, . . ., Qt corresponding to x\, . . . ,x n , the probability of obtaining 
f(x\, . . . , x n ) when measuring the first qubit of 

UtQtUt-1 ■ ■ ■ UiQ T U \ll> start) 

is at least 2/3. 
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